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Abstract: We minimize the one-loop effective potential for SU{N) gauge theories in- 
cluding fermions with finite mass in the fundamental (F), adjoint (Adj), symmetric (S), 
and antisymmetric (AS) representations. We calculate the phase diagram on x as a 
function of the length of the compact dimension, /3, and the fermion mass, m, for various 
N and A'^^. We consider the effect of periodic boundary conditions [PBC(-|-)] on fermions 
as well as antiperiodic boundary conditions [ABC(-)]. With standard ABC(-) on fermions 
only the deconfined phase is found at one-loop for all representations considered. However, 
the use of PBC(-I-) produces a rich phase structure. These phases are distinguished by the 
eigenvalues of the Polyakov loop P. In the case of fundamental representation fermions 
[QCD(F,+)], a phase in which Re TrP is minimized (and negative) is favoured for all values 
of mf3. For N odd charge conjugation (C) symmetry is spontaneously broken in this phase 
due to 0{1/N) effects. Minimization of the effective potential for QCD(AS/S,+) results 
in a phase where |ImTrP| is maximized, resulting in C-breaking for all N and all values 
of m/3, however, the partition function is the same up to 0{1/N) corrections as when 
ABC are applied. Therefore, regarding orientifold planar equivalence, we argue that in 
the one-loop approximation C-breaking in QCD(AS/S,+) resulting from the application of 
PBC on fermions does not invalidate the large N equivalence with QCD(Adj,-). Similarly, 
with respect to orbifold planar equivalence, breaking of Z2 interchange symmetry resulting 
from application of PBC to bifundamental (BF) representation fermions does not invali- 
date equivalence with QCD(Adj,-) in the one- loop perturbative limit because the partition 
functions of QCD(BF,-) and QCD(BF,+) are the same. Of particular interest as well is 
the case of adjoint fermions where for Nf > 1 Majorana flavour confinement is obtained 
for sufficiently small m/3, and deconfinement for sufficiently large mf3. For A > 3 these two 
phases are separated by one or more additional phases, some of which can be characterized 
as partially-confining phases. 
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1. Introduction 

There has been substantial recent interest in properties of QCD-like gauge theories. These 
field theory models typically involve SU(N) gauge theories with fermions in representations 
beyond the fundamental, including the adjoint, symmetric, and antisymmetric representa- 
tions. In addition to the standard antiperiodic boundary conditions applied to fermions 
there are also good reasons for studying cases in which periodic boundary conditions are 
applied. For example, SU{N) gauge theories with adjoint fermions, and related theo- 
ries that preserve the center symmetry and confining properties of the pure gauge theory 
||l|, ^, have proven to be very useful in understanding mechanisms of confinement 

lH> 05 §5 0, 0]- In particular, lattice simulations and analytical results both indicate the 
existence of a perturbatively accessible confining region, which is analytically connected to 
the confined phase of the pure SU{N) gauge theory, when periodic boundary conditions 
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are applied to an adjoint potential contribution |2|, |T^. These theories are also useful in 



the study of Eguchi-Kawai reduction |13], in which confined gauge theories should exhibit 
volume independence in the large N limit Q. 

Gauge theories with fermions in higher dimensional representations than the fundamen- 
tal, such as the two index symmetric, antisymmetric, and adjoint representations are also 
required in the equivalence between a supersymmetric theory and a non-supersymmetric 
one, in the large N limit, known as orientifold planar equivalence |14, These models 



are also relevant for the construction of conformal ITq, and near-conformal |18, 19 1 



QCD-like gauge theories |20|, ^, |2^, |23|, gj, |25|, g6||27|, g8|. 



It has been known for some time that perturbation theory is not valid for calculation 
of phase transitions from a weakly-coupled phase to a strongly-coupled phase which is clear 
from high temperature calculations in pure gauge theories and QCD |2^. However, various 
phases are accessible in QCD-like theories at weak coupling when PBC(-I-) are applied to 
fermions. For example, when fermions are in the adjoint representation perturbation theory 
has revealed an exotic phase structure, and C-breaking phases are observed for fundamental, 
symmetric or antisymmetric representation fermions. 

It is possible to calculate the phase diagram using high temperature perturbation the- 
ory for SU{N) gauge theories with fermions in an arbitrary representation. We refer to 
these theories as QCD(R), where R is the representation of fermions used. QCD(F,-) cor- 
responds to ordinary QCD when = 3, where the (-) indicates antiperiodic boundary 
conditions (ABC), in the time dimension, on fermions. When considering periodic bound- 
ary conditions (PBC) on fermions we use the symbol The gauge fields will always 
have the usual periodic boundary conditions applied. Our results reside on the topology 
X where the time dimension is compactified with length (3. When ABC are applied 
to fermions, this corresponds to taking the theory at finite temperature. ^ 

In this paper we provide a catalog of some of the phases observable from one-loop 
perturbation theory in QCD(R,ib) for a single small compact dimension with fermions in the 
fundamental (F), symmetric (S), anti-symmetric (AS), and adjoint (Adj) representations. 
The one-loop calculations of In Zpi^± for finite fermion mass m, finite chemical potential 
and constant finite gauge field A4 are included in Appendix ^ where we provide a detailed 
derivation which compiles results scattered among various sources: For pure SU{N) gauge 
theory, V^fj has been calculated to one- loop order in [^. In |^, ^ the contribution 



of fermion mass is included at one- loop for finite T and ^. In |32] the one and two- loop 
results are calculated for QCD with one flavour of massless fermions and finite /x. In |36| 
the one-loop effective potential was calculated in QCD(AS/A/Adj) for the case m = and 
/X = 0. 

The partition function at weak coupling has also been clearly derived at one-loop for 
SU{N) gauge theories on x in the very clear and pedagogical papers [33, In |33| 
the authors calculate the large-A^ deconflning phase transition temperature for the pure 



^When PBC are applied to fermions the theory does not correspond to a finite temperature field theory 
[ABC are required in that case due to the trace over (physical) anticommuting fields in Z — Tr (e'''^)]. 
However, the formalism is the same so we do not refrain from using the language of thermodynamics, even 
though it only strictly applies in the case of ABC on fermions. 
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Yang- Mills theory (and the Hagedorn transition temperature in = 4 SYM) on x S^. 
In 1 34] they show that the Yang- Mills theory deconfinement transition is first-order. Gauge 
theories on x and x have different, although slightly overlapping regions of 
perturbative validity so it is important to study both. To obtain true phase transitions on 
a finite manifold like x S^, it is necessary to take the large N limit. On x M^, this 
is not necessary, but perturbative calculations on x with a small number of fermion 
flavours ^ are only valid in the limit of small where Rgi <C ^qcd (see also the end of 
1 34] for a nice discussion of this point). 



2. Conventions and methods 

All of the thermodynamics of QCD(R) can be obtained from the partition function ^QCD(i?)- 
For Nf quark flavours in representation R, with masses mj, at chemical potential and 
at inverse temperature P = 1/T, the QCD(R) partition function is given by the (Euclidean 
space) path integral 



where the gauge field A^{x) = T^A^(2;), ^ = 0, 1, 2, 3, the a indices are a = 0, A^^ — 1, 
and the are the generators of SU{N) in the representation R. ^ is a D{R) component 
vector of fermion fields in the representation R containing Nf anti-commuting 4-spinors 
ipf. ip contains the corresponding antifermion fields. M is the fermion mass matrix where 
{M)ff' = rafSffi and M is the fermion chemical potential matrix where {M)ffi = ^ifSff. 
~ Tm-^m where Dfj^ is the covariant derivative 

D^, = d,, + A^. (2.2) 

The are D{R) x D{R) antihermitian matrices that transform as the representation R 
of the fermion fields on which they act. 

Since our interest is in the phase diagram of this theory we introduce the Polyakov 
loop order parameter which is defined as the path-ordered exponential of the temporal 
component of the gauge field, 

P(f) =Pe^o''dt^oW_ (2.3) 
For a constant background field defined hy Aq = iv/(3 the Polyakov loop is 

p = e^^o = (2.4) 

where for simplicity we have chosen a gauge in which Aq is diagonal and v is real, diagonal 
and traceless with elements {v)ij = Vi5ij. In this case we take Aq (and thus v), which 
transforms as the representation R, to be in the form of an x matrix. Then 

^For Nf — 5 Majorana flavours the validity of calculations on x R"^ extends also to the limit of large 
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P = diag{e™S...,e™^}. 



(2.5) 



2.1 Effective potential 

Details of the derivation of In Z on S*^ x for arbitrary fermion mass m and chemical 
potential fi is included in Appendix Using this result we obtain the one-loop effective 
potential in terms of the Polyakov loop, 



Vi-ioopiP,m,f3,n)± 

■lnZ(P,m,/3,/i)± 



I3V3 



1 



PV3 



[-2A^/lndet (-D|(P) + m^) + Indet {-D\aj{P))] 



2 AT 



7r2/32 



■E 

n=l 



(±1)^ 



2 A 1 



en/3/.Tr^(ptn) + e-"^^Tr^(P") K^infim) - ^ V —Tta{P 



n=l 



+ 



1 ^ 2 



(2.6) 

where TVyi indicates a trace over an object in the adjoint representation, and gR,± depends 
on the group representation of the fermions and the boundary conditions. It is defined as 
the trace of the matrix 



inv—n(3(mt— fj) 

'^=1 (2.7) 



aP(mt-n) _ 2 cos V + e^^(™*^^) 



in the representation R. The trace in various representations is determined using the Frobe- 
nius formula combined with tensor product methods as discussed in Appendix ^ We define 



= TrR g{l3, mt, v) (2.8) 

for periodic boundary conditions applied to fermions. The values of gR^± for both PBC 
and ABC applied to fermions can be found in Appendix ga^- are obtained from gii^+ 
by taking v ^ v + tt. 

In what follows we consider the phase diagram of QCD(R) for fermions in the Fun- 
damental (F), symmetric (S), antisymmetric (AS), and adjoint (Adj) representations. We 
take = so the result simplifies to 
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_ 2m'^Nf ^ (±1)" 



7r2/32 

n=l 



(+1)" 2 1 

^ [Tr^P")] K,{nPm) - ^ ^ ;^Tr^(P^ 



n=l 



dt(t2 - l)3/2Re [g^^^ (/5,mt,0,r;)] 



1 



247r2 



AT 

El 

«J=1 



(27r 



IT 



1) 



2.2 Chiral condensate 



The chiral (or quark) condensate, {ipip), is an order parameter for the chiral symmetry of 
a theory. limm^o(V'V') 7^ indicates that chiral symmetry is broken. It also serves as an 
indicator of the order of transitions between phases. From the effective potential it is easy 
to get the chiral condenstate, which is given by the mass derivative 



{tpil))i-loop{'m) 



lim ^ \nZ{m) = ^-^Veff{P,m). (2.10) 

Vi^oo V^Nf dm ^ ' Nf dm ' ^ ' 



With this definition in QCD(R) for fixed (''/'V') is independent of N ^ if the effective 
potential is minimized in the same phase (the fermion contribution to V^// is proportional 
to Nj^ and this is the only term with mass dependence). In the following sections we 
measure (V'V') along with the effective potential. 



3. Results 

To see clearly the effect of fermion mass m on the phase diagram of QCD(R) we perform 
calculations at ^ = 0. In this section we numerically minimize V^jf in eq. ( |2.9[ ) for m/? from 
to 10. To ensure accuracy of the results we use a smaller stepsize A(m/3) = 0.1 in the 
region < mf3 < 3 since this is where we observe the most phase transitions. The use of a 
smaller step size also serves as a check that the minimization routine is finding the global 
minimum. We minimize 14// using an algorithm that searches for the minimum using 10 - 
200 random starting points (where more were needed for larger and complicated phase 
diagrams) to increase the probability that the global minimum is obtained rather than just 
a local minimum. 



3.1 Fundamental fermions 

For SU{N) gauge theories with fundamental representation fermions, to which either an- 
tiperiodic (-) or periodic (+) boundary conditions have been applied, there are three pos- 
sible phases which are distinguished in Table |^ according to their Polyakov loop eigenvalue 
angles v = {vi, ...,vj\f}, as well as TvpP. 

The confined and deconfined phases are familiar. The confined phase is that of pure 
gauge theory. It is located in a region of strong coupling and is thus only accessible non- 
perturbatively, for example by using lattice simulations. QCD(F) corresponds to pure 
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Table 1: Table of the possible phases of SU{N) gauge theory with fundamental representation 
fermions to which either antiperiodic (-) or periodic (+) boundary conditions have been applied, 
distinguished according to their Polyakov loop eigenvalue angles v = {vi, ...,vn}, as well as TipP. 
Note [n, m] = n mod m. 



Phase 


V 


TrpP 


confined 

deconfined 

anti-deconfined 


riV-l N-3 N-^ttI 
{0,0, ...,0} 

{ TT ± f [iV, 2] , vr ± f [iV, 2] , . . . , TT ± f [iV, 2] } 




N 

-iVexp(±if [iV,2]) 



SU{N) gauge theory when the fundamental fermion mass m oo. In the case of funda- 
mental fermions with non-infinite mass at a value of /5 = 2N/g'^ for which the confined 
phase is observed in the pure gauge theory, we expect v and TipP to be different from 
those of the confined phase in Table || ^. This is because fundamental fermions explicitly 
break the Z{N) center symmetry of the theory. However, v and TrpP should approach 
the pure gauge values as m — > oo. 

3.1.1 Antiperiodic boundary conditions [ABC(-)] 

Minimizing the one-loop effective potential of eq. ( |2.9| ) we find that when ABC are applied 
to fundamental fermions the deconfined phase of Table |l| is always favoured. To consider 
a physical theory, we take = 3 in addition to ABC on fermions, which gives QCD at 
finite temperature. As expected in the perturbative limit, for all m(3 the effective potential 
is only minimimized for Polyakov loop angles corresponding to the deconfined phase, such 
that TipP is magnetized along the positive real axis: 

v = {0,0,0}; TrirP = 3. (3.1) 

The effective potential is plotted in Figure || (Left). In this type of plot the dots 
correspond to the minimization of eq. ( |2.9| ) with respect to the Polyakov loop angles Vi for 
a range of m/3 from to 10. The curves correspond to possible phases of QCD distinguished 
by the values of the Polyakov loop angles as presented in Table |l[ Even though only the 
deconfined phase is accessible via perturbation theory in QCD, a significant feature in 
Figure |^ (Left) is the presence of the inflection point in 14// at mfi « 1.4. This implies a 
large one-loop contribution to the chiral condensate (V'^) as indicated in Figure |l] (Right). 
Including non-perturbative contributions (^V') = is expected in the high temperature 
limit of QCD. 

3.1.2 Periodic boundary conditions [PBC(-|-)] 

When PBC are applied to fundamental fermions Vp(/^£)(_|_) always favours the anti-deconfined 
phase of Table |l] as indicated in Figure ^. In this case the Polyakov loop angles take the 

^This behaviour is suggested in js^ ] using 1-loop results for SU{N) gauge theory with both adjoint and 
fundamental representation fermions (the adjoint fermions are used to make the confined phase accessible 
perturbatively) . The confined phase has v = {0, (j), —(j}}. When the fundamental fermion mass mp — » oo, 
4> — > 27r/3 as expected. As mp is decreased from infinity (j) also decreases. 
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Figure 1: QCD: (Left) Vfund{-) for ^/ = 1; (Right) (V-V) FUNDI-)- 





Figure 2: Vpijis[u(+) for Nf — 1. The dots correspond to minimization of the effective potential 
of QCD(F,+) with respect to the Polyakov loop angles Vi. The cm^ves correspond to evaluation of 
the effective potential for Vi = tt ± [N, 2] , indicating that the anti-deconfined phase of Table |l| 
is favoured for all m(3. The curves corresponding to other phases in Table |l] are not included for 
clarity, and because they would correspond to larger values of the effective potential. 



values which are closest to v = {vr, vr, vr}. As indicated in Table 0, for even the 
Polyakov loop angles take those exact values such that Tr^^P = —N, but for N odd they 
can only come close, taking the values 

Vi = TT±^ Vi, (3.2) 

which makes P complex. In other words, for N odd the Polyakov loop eigenvalues will take 
the value of either of the two Nth roots of unity which are closest to — 1 . In this case the C- 
symmetry is broken since the Polyakov loop is not invariant under P ^ P* . For the familiar 
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case of = 3 but with PBC on fermions either of the two complex phases corresponding 
to Vi = ib27r/3 are preferred, whereas for ABC the phase with = is preferred. Ah three 
correspond to deconfined phases of the pure gauge theory. Therefore, when fundamental 
fermions are added to pure gauge theory in a deconfined phase, preference of a particular 
vacua depends on the type of boundary conditions applied to the fermions. 



3.2 Antisymmetric (AS) / Symmetric (S) fermions 

For SU{N) gauge theories with fermions in the antisymmetric (AS) or symmetric (S) 
representation, the possible phases are similar. In the non-perturbative regime, for ABC 
on fermions, we expect to find a phase which approaches the confined phase in the m — > oo 
limit. This phase is not observed in our one-loop perturbation theory calculations. As 
expected, for AS/S fermions with ABC, the deconfined phase with v = {0,0, ...,0} always 
minimizes Veff- With even the deconfined phase can also be defined by v = {vr, vr, vr} 
since the effective potential only includes terms where the eigenvalues are subtracted (boson 
contribution) or added (fermion contribution). However, with PBC on fermions a phase in 
which the Vi are all the same, and as close as possible to ib7r/2 is favoured. This is because 
for tensor representation fermions, gAS{+) ^'^d gs{+) of eqs. ( A.2g| ) and ( A.29|) depend on 
Vi + Vj, where in the case of fundamental fermions gF{+) of eq. ( A.26|) depends on 



Vi 



where values closest to vr were preferred. Since Vi close to ib7r/2 is preferred for all when 
considering AS/S fermions with PBC, the C-symmetry is always broken [the exception is 
N = 2 QCD(S) as it is equivalent to QCD(Adj)]. 

For = 4, 8, 12, ... the C-breaking phase has 2 possible vacua corresponding to 

vr 

Vi = ±- Vi. (3.3) 

Notice that these are exchanged under C. 

For A^ odd the (double) C-breaking phase has 4 minima, 

= =^ (i * 27?) P'"' 

where '^, = ±(^ + 2^) are exchanged under C, and Vi = ± — are exchanged under 
C. 

For A^ = 2, 6, 10, ... the (double) C-breaking phase has 4 minima. 



t;, = ± (- ± -) Vi, (3.5) 

where Vi = ±(^ib-^) are exchanged under C, and fj = ±(| — ■^) are exchanged under 
C. The effective potential for various A^ minimized in the C-breaking phases is shown in 
Figure | (Left) for QCD(AS,+) and Figure | (Right) for QCD(S,+). The possible phases 
of SU{N) gauge theory with AS/S representation fermions and either ABC or PBC on 
fermions are summarized in Table 
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Table 2: Table of the possible phases of SU{N) gauge theory with antisymmetric (AS) or symmetric 
(S) representation fermions and either antiperiodic (-) or periodic (+) boundary conditions on 
fermions. [n, m] = n mod m. 



Phase 


V 


confined 

deconfined 

C-breaking 


fN-l N~3 N-1 ^ 
1 TV ' AT 1 N i 

{0,0, ...,0} (for N even this phase can also be defined by {'/r,7r, ...,7r}) 
V, = ±[N,2] (i ± 2^) ± i[iV,4][(iV + 1),2] (f ± f ) 
±i[(7V + 2),4][(iV + l),2](f) 




Figure 3: Nf = 2: (Left) Vas{+)', (Right) Vsymm{+)- The dots correspond to minimization 
of Veff in eq. (2^) with respect to the Vi. The curves result from evaluation of eq. (2^) for 
V, = ±[N,2] (f ± ± i[7V,4][(iV + 1),2] (f ± f ) ± i[(7V + 2),4][(7V + 1),2] (f), as in Table |. 
This indicates that in QCD(AS/S,+) for any the C-breaking phase is favoured for all m/3, except 
for A = 2 QCD(S), which is equivalent to QCD(Adj) [see also Figure §. 



3.3 Adjoint fermions 

QCD(Adj) is unique among the models studied here in several regards. For one, quantities 
including only adjoint representation traces over the Polyakov loop, Tr^^P, like the partition 
function ZAdj, ai^e clearly invariant under Z{N) transformations [Tr^^P = |Tri?P|^ — 1 is 
invariant under P — > zP for z £ Z{N)]. This means that all vacua which are Z{N) rotations 
of each other are equivalent, as is the case for the pure gauge theory. Since all phases have 
at least one of the vacua lying on the real axis, we never observe C-symmetry breaking 
in QCD(Adj) for either PBC or ABC applied to adjoint fermions. In addition, the phase 
diagram of QCD(Adj,-|-) with PBC on fermions becomes quite rich for Nf > 2 Majorana 
flavours {Nf Majorana flavours = 2Nf Dirac flavours). For QCD(Adj) we take Nf to be 
the number of Majorana flavours since adjoint fermions are their own antiparticles. This 
means that Nf Nf/2 in eq. (p|). 

When ABC are applied to adjoint fermions the deconfined phases are always favoured 
in the minimization ofVf,ff. In the deconfined phases the Vi are all the same and correspond 
to one of the Nth roots of unity, as is the case for the pure SU{N) gauge theory. 

Changing the fermion boundary conditions to periodic we find that for Nf = 1 the 
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Figure 4: (Left) Vadj(+) with Nf — 1 Majorana flavour. The dots correspond to minimization of 
Veff in eq. (2.9) with respect to the Vi. The curves correspond to the result of evaluating eq. (2.9) 
for Vi = 0V«, the values corresponding to the deconfined phase. (Right) {^^) adj{+) for Nf = 1 
and various values of N. 



deconfined phases are also favoured for all N as shown in Figure ^ (Left). The correspond- 
ing calculation of the chiral condensate, {^^)i^ioop for QCD(Adj,+) is shown in Figure ^ 
(Right). However, for Nf > 2 there is a great variety of observed phases. In particular the 
confined phase appears to be one of the phases present for all N. The observed phases of 
QCD(Adj,+) for Nj >2 and from 2 through 9 are given in Tables | - |l^. Recall that 
Z{N) rotations of Polyakov loop eigenvalues included in the tables result in the same value 
for Vf.ff since it is Z(A^)-invariant for QCD(Adj). 
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Table 3: Table of the observed phases N = 2 QCD(Adj,+) for Nf >2 Majorana flavours. 



Phase 


V 


TrpP 


confined 
deconfined 


r TT TT \ 

12' 2 J 

{0,0} 



2 



Table 4: Table of the observed phases N = 3 QCD(Adj,+) for Nf > 2. 



Phase 


V 


TvpP 


confined 


I 3 ' ^' 3 J 





split 


{0,7r,7r} 


-1 


deconfined 


{0,0,0} 


3 



Table 5: Table of the observed phases N = 4 QCD(Adj,+) for Nf > 2. 



Phase 


V 






confined 


f 3n TT TT 3tt\ 
I 4 1 4 ' 4 ' 4 J 








S'?7(2)-confincd 


J TT TT TT TT\ 

I 2 ' 2 ' 2 ' 2 J 





-4 


deconfined 


{0,0,0,0} 


4 


4 



Table 6: Table of the observed phases iV = 5 QCD(Adj,+) for Nf > 2. 



Phase 


V 




confined 

attractive 

split 

deconfined 


/ 47r 27r PI 2tt Att \ 
t 5 ' 5 ' ^' 5 ' 5 J 

{0,-0,-0, 0, 0} 

{0,0,0,7r,7r} 
{0,0,0,0,0} 




1 + 4cos(0) 

1 

5 



Table 7: Table of the observed phases iV = 6 QCD(Adj,+) for Nj > 2. 



Phase 


V 






TyfP^ 


confined 


/ 57r Stt TT TT 37r Stt \ 
16'6'6' 6' 6' 6-1 











SU (3)-confined 


/ 27r PI 2tt 2tt n. 2tt \ 
I 3 ' ' 3 ' 3 ' ' 3 J 








6 


5C/(2)-confined 


t TT TT TT TT TT 7r\ 

l2' 2'2' 2'2' 2J 





-6 





deconfined 


{0,0,0,0,0,0} 


6 


6 


6 



Table 8: Table of the observed phases N = 7 QCD(Adj,+) for Nf > 2. 



Phase 


V 


Tr^P 


confined 
repulsive 
? 

attractive 
split 

deconfined 


/ Gtt 4-7r 2tt r\ 2tt Att Qtt \ 
17)777)'-') 7) 7) 7/ 

{0,-0,-0,0, 0,7r,7r} 
{0,0,0,-0, -0,0,0} 

{0,-0,-0,-0,0,0,0} 
{0, 0, 0, TT, TT, TT, tt} 

{0,0,0,0,0,0,0} 




-1 + 4cos(0) 
3 + 4cos(0) 
1 + 6cos(0) 
-1 
7 
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Table 9: Table of the observed phases iV = 8 QCD(Adj,+) for Nf > 2. Here C(0) = cos(6'). 



Phase 


V 


TrpP 


TrirP2 




confined 


r Ttt StT StT TT TT Stt Stt Ttt t 
t8'8'8'8' 8' 8' 8' 8J 











Sf/ (4) -confined 


r Stt tt tt 3tt Stt tt tt 3tt \ 
1454' 4' 4'4'4) 45 4/ 








-8 


attractive 


{0,0,-(p,-(j),-(p,(j),(f),(f)} 


2 + 6c(^) 


2 + 6c(2<^) 


2 + 6c(4^) 


5C/(2)-confined 


f TT TT TT TT TT TT IT TT 1 

l2' 2'2' 2'2' 2'2' 2J 





-8 


8 


deconfined 


{0,0,0,0,0,0,0,0} 


8 


8 


8 



Table 10: Table of the observed phases N = 9 QCD(Adj,+) for Nf > 2. Here C(e) = cos{e). 



Phase 


V 




Tr^P^ 


confined 
mixed 

SU (3)-confined 

attractive 

split 

deconfined 


/ Stt 67r ri 67r 87r \ 
t 9 ' 9 ' 9 ' 9-1 

{0, -(/>, -(j), (p, (f), -X, -X, X, x} 

r 27r PI 47r 2tt r, Itt 2tt n Att \ 
13'' 3 ' 3 ' ' 3 ' 3 ' ' 3 J 

{0, -(f), -0, -(f), (f), (f), (f), (f)} 

{0, 0, 0, 0, 0, vr, TT, TT, tt} 
{0,0,0,0,0,0,0,0} 




1 + 4[C(^) + C(;,)] 



1 + 8cos(0) 

1 

9 




1 + 4[C(30) + C(3x)] 

9 

1 + 8cos(3^) 

9 

9 
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The results from minimization of the effective potential for QCD(Adj,+) with respect 
to the Polyakov loop angles Vi for Nj = 2 and N from 2-9 are shown in Figures ^ - 12 



These phase diagrams are quite rich and appear to increase in complexity as A*" increases. 
The confined phase is defined by 

v = {^^^,^^vr,..., ^vr}; Tr^P = 0. (3.6) 

The deconfined phases are defined by 

v = {0,0,...,0}; TtfP = N, (3.7) 

and Z{N) rotations. 

For N not prime there are additional phases which correspond to confinement of a 
subgroup L of SU{N), which are called S'C/(L)-confined In these phases quarks are 
still confined, however groups of L quarks are not. The Polyakov loop angles correspond 
to N/L copies of the L angles of the SU{L) confined phase. For example, in 5C/(6) we 
observe SU{2) and SU{3) confined phases. The SU{2) confined phase corresponds to 



vrTrTrvrvrTT 
^2'~2' 2'~2' 2'~2^' ^^'^ 



and the SU (3) confined phase has 

,27r 27r 27r 27r, , , 

as indicated in Table 0. These phases can be distinguished by taking the trace of powers of 
the Polyakov loop. In the SU{2) confined phase Tr^^i-" = TrpP^ = 0, however TrpP'^ ^ 0. 
In the 5;7(3) confined phase Tr^P = Tr^P^ = q, but Tr^P^ ^ 0. The same pattern holds 
for arbitrary L. 

In addition, for = 5,7,8,9 we find an attractive phase, called such because TrpP 
is not constant with respect to small changes in m/?, rather ITtjt'PI slowly decreases as 
mP increases, indicating that the Polyakov loop eigenvalues are attracted together. In the 
repulsive phase of SU{7) the Polyakov loop eigenvalues are repelled apart with increasing 
m(3. In SU(9) there is a mixed phase where some of the eigenvalues attract, while others 
repel. 

For N odd we also find split phases which favour a small magnitude for |Tri?P|. These 
have {N it l)/2 angles that are vr and {N =f l)/2 angles which are 0, where the top sign 
corresponds to N = 3,7,11,... and the bottom sign corresponds to N = 5,9,13,..., such 
that det P = 1 as required. The Z{N) symmetry is completely broken in these phases, 
albeit in a special way. For example, in SU{5) this phase is given by v = {0, 0, 0, vr, vr} 
which is an SU{3) x SU{2) [mod Z{5)] phase. 

In Q we observed the N = 3 split phase in lattice simulations of a related theory and 
referred to the phase as "skewed" in that Pro] z(n)['^^fP] < as opposed to the deconfined 
phase where Proj2(Ar) [Tr^^P] > 0, causing the vacua of the skewed phases to lie at angles 



* Partially confined phases are also discussed in ^ 
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midway between those of the deconfined phases. However, for = 5, 9, 13, ... the vacua of 
the spht phases have angles corresponding to the N roots of unity, as do the deconfined 
phases. As will become important shortly, the split phase of odd N and the S'C/(2)-confined 
phase of even are related in that Vj^dj^^'j is the same in these phases, in terms of N. 

It is important to note that additional phases are also possible. Our step size in the 
region 0.1 < m-/3 < 3.0 is A(m/3) = 0.1, and for 3 < m(3 < 10 the step size is A{m(3) = 1. 
So there is room for additional phases in narrow regions of m/3 that extend less than 0.1 
for mf5 < 3.0, and less than 1.0 for 3.0 < m/? < 10. 



3.3.1 Large N limit 

From the phase diagrams of QCD(Adj,+) in Figures ^ - ^ it is clear that the confined 
phase extends for a smaller range of m(3 as increases. To determine if the confined phase 
persists into the limit of very large we minimized in the region of small m/3 for 
values of A^ up to 19. The value of m/3 at which there is a transition out of the confined 
phase was measured. The results are plotted in Figure 13. The curve represents a lower 
bound: for each value of A^ on the curve the value of m/? 0.1 above it does not belong 
to the confined phase. The confined phase becomes quickly less accessible as A^ increases 
from 2 to 10, but appears to level off in the large A^ limit providing a narrow region, near 
m/3 ~ 0, in which it is the preferred phase. 

Another noteworthy feature of the phase diagrams in Figures |5| - 12 is that the tran- 
sition to the deconfined phase always occurs at m/3 = 2.00199, independent of A^. This is 
explained in the next subsection. 



3.3.2 Variations in Nf 

We also calculated the phase diagrams in QCD(Adj,-|-) for Nf = 3 flavours of adjoint 
fermions for A'^ = 3,4,5,6. The general trend is that the value of m/3 below which the 
confined phase is preferred, given by {m(3)criti increases with Nf. This is illustrated in 
Figure 0. 

Overall, increasing A'^^ can offset the decrease in [mf3)crit with A^, but there is a natural 
limit in the value of Nf in that to preserve asymptotic freedom in QCD(Adj) it is necessary 
to keep Nf < 5 Major ana flavours. 

Also, as A'^^ increases the transition to the deconfined phase, which is still independent 



of A^, also increases, as shown in Figure 15. The value of mf3 at which this transition is 
realized is given by the solution of 



E 



1 



{ (2n 



NfK2 [(2n - l)m(3] 



[{2n-l)m(3Y 



0. 



(3.10) 



As mP is increased the transition to the deconfined phase always occurs from an SU{2)- 
confined phase (for A^ even) or a split phase (for A^ odd). The formulas for VAdj{+) 
identical in the split and 5'C/(2)-confined phases since in both cases Tr^(P") = for n 
odd and TrA(P'') = N^ - 1 for n even (see Tables | - |l0D. Since TrA(P") = A^^ _ ^ 
for all n in the deconfined phase the only terms that differ when comparing the split, or 
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S'C/(2)-confined phase, and the deconfined phase, occur for n odd. Setting the formulas for 
VAdj{+) ill the spht (or 5J7(2)-confined) phase and the deconfined phase equal then gives 
the equality in eq. ( 3.10| ). 



3.3.3 Chiral condensate 

The chiral condensate in QCD(Adj) is given by 

where we replaced V'V' in eq. ( |2.10| ) with A A since adjoint particles and antiparticles are 
indistinguishable. In order to ascertain the order of the transitions in our phase diagrams 
of VAiij{+) vs. m(3 in Figures ^ - 12, we calculated also the chiral condensates. These 



are shown for N from 2 through 9 in Figures - The most remarkable transition in 
all cases is that to the deconfined phase which exhibits the largest jump in m/3. But, all 
phase transitions are marked by a noticeable jump in (AA) suggesting that all transitions 
are first order. There is also a smooth increase in (AA) at low m/3 as observed in the 
case of fundamental fermions, and as well for symmetric and antisymmetric representation 
fermions. 

3.4 Orientifold and orbifold planar equivalence 

In 2003 - 2004, Armoni, Shifman, and Veneziano published a series of papers in which they 
proved non-perturbatively the equivalence of the bosonic sectors of QCD(Adj) with Nf 
Major ana fermion flavours and QCD(AS/S) with Nf = Nf irac fermion flavours, in the 
planar (or large N) limit This is called orientifold planar equivalence. For = 3, 

QCD(AS) is equivalent to QCD (taking ^ Also, QCD(Adj) with one Majorana 

flavour of massless fermions is = 1 supersymmetric Yang-Mills (SYM) theory. This 
makes orientifold planar equivalence particularly useful since it allows one to use results in 



supersymmetry to make estimates in one-flavour QCD up to 0{1/N) corrections [38]. For 



example, in [39| the quark condensate was calculated in one-flavour QCD from the gluino 
condensate in SYM. This was subsequently verified in using lattice simulations. 

Perturbative comparison of the phase diagrams in QCD(Adj) and QCD(AS/S) was per- 
formed first in [^], then [^ ], then quite thoroughly in |^2|. In |^6| it is argued that orien- 
tifold planar equivalence is only valid when charge conjugation (C) symmetry is not broken. 
In the example provided, they showed that C-symmetry is broken in U{N) QCD(AS/S) 
with periodic boundary conditions for massless fermions on x M^, for small S^. While 
it seems to be true that C-symmetry is required in general for planar equivalence, their 
example is a special case in which orientifold planar equivalence actually still holds, despite 
C-parity being broken. From the derivation of the fermion determinant in Appendix it 
is clear that changing fermion boundary conditions from periodic to antiperiodic results in 
shifting the gauge field by iir/ (3, i.e. v — > v + tt. For the U{N) theory (which corresponds to 
the SU{N) theory for A^ = 4, 8, 12, ...) transforming from antiperiodic to periodic bound- 
ary conditions results in the the effective potential being minimized when the individual 
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Polyakov loop angles are shifted each according to Vi ^ Vi + TT/2. In terms of the partition 
function, 



QCD{AS/S,+) : lnZAs/s,+ oi^i + vj), Vi = 7r/2, 

QCD{AS/S, -) : In Z^s/s,- ~ 9{vi + f j + vr), v, = f). ^ ' ' 

Therefore, for the U{N) QCD(AS/S) theories the partition function will always be the same 
under transformation of fermion boundary conditions between periodic and antiperiodic. 

In summary the equivalence holds in the perturbative limit when considering QCD(Adj,- 
), QCD(AS,+), and QCD(AS,-), but not QCD(Adj,+). This is shown explicitly in Figure 
(Left) where the one-loop effective potentials of these theories in U (N) are plotted to- 
gether for N = 8. QCD(AS,+) and QCD(AS,-) correspond exactly. QCD(AS,+/-) and 
ACD(Adj,-) correspond up to 0{1/N) corrections. These statements of course also hold 
for QCD(S) in place of QCD(AS). For the SU{N) theories, the partition functions for 
QCD(AS/S) agree exactly with the U{N) case when is a multiple of 4, otherwise they 



agree up to (D{l/N) corrections, as seen in eqs. (gj), (U). In QCD(Adj), the SU{N) 



theory has one less degree of freedom corresponding to an extra —1 in Tr^P = iTrjT'Pl^ — 1 
(see Appendix^, resulting in an extra constant in gAdj,± in the partition function [see eqs. 
( |A.27] ) and (A.31)]. This constant is subleading in the large N limit. Therefore, the above 



statements about orientifold planar equivalence hold in both the U{N) and the SU{N) 
theory, as expected. 

Conclusions analogous to those for orientifold planar equivalence can also be shown 
to hold in a related orbifold planar equivalence as well. In |43| Kovtun, Unsal and Yaffe 



proved non-perturbatively the large-N orbifold equivalence of [U{2N)] QCD(Adj) with N^^ 
Majorana fermion flavours and QCD(BF) [U{N)i x U{N)2 gauge theory with bifundamen- 
tal representation fermions] with NP = N¥ Dirac fermion flavours, in the large mass. 



strong coupling phase. This result was anticipated in |44, The technique of orbifold 



projection is quite useful and has been of recent interest as a means of testing volume 
independence in confining gauge theories via Eguchi-Kawai reduction [46|. 



In [47, 48, ^ the breaking (or lack thereof) of Z2 interchange symmetry of the two 



gauge fields is explored, the symmetry which is thought to be required for orbifold planar 



equivalence to hold. In [47| this Z(2)-symmetry was shown to be broken in QCD(BF,-|-) 
with a compact dimension. In [^] it was argued that this symmetry is broken on 
when periodic boundary conditions are applied to bifundamental representation fermions, 
causing the equivalence to fail. In [49| it was argued that there is no evidence to support 



Z(2)-breaking on R^. 

As with the orientifold case, it is straightforward to show perturbatively that the effect 
of symmetry breaking resulting from the application of PBC on bifundamental fermions 
does not change the partition functions from that of the theory with antiperiodic boundary 
conditions on fermions ^. The effect of Z(2)-symmetry breaking in QCD(BF,-|-) is cancelled 
out by the application of PBC on fermions so the equivalence does not fail. But, if the 



^We would like to thank Mithat Unsal for pointing out this important calculation to us. 
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Z(2) symmetry is broken by some means other than apphcation of PBC on fermions, such 
that the partition function is different, then the equivalence would not hold. 

QCD(BF) is a U{N)i x U{N)2 theory with gauge fields in the adjoint representation of 
each group, and bifundamental representation fermions which transform as the fundamental 
in the first group, and antifundamental in the second. The additional gauge field results 
in a doubling of the gauge boson term, which results from the action 



QCD{BF) 



1 

4^ 



dr / d-'xTiF {Fi,^yFi^^u + F2,f,^F2,f,u) + {iPsF " loM + M) ^. 

(3.13) 

D^p = + with A^^ = A^ + A2 where Ai transforms as the fundamental and A2 
transforms as the antifundamental so the Polyakov loop is P = e*^ with v = vi + V2. The 
fermion contribution to the effective potential follows from eq. ( [A. 24 ) for the determinant 
with fermions in representation R = BF so TtbfP = Tri;'PiTri;'P2 S'^d Pj = e^'"^ = 
(iiag{e™^'^ , e*''^-'^}, j = 1,2. The effective potential (for = 0) has the form 



Vbf,-. 



1 



2m^Nf 
37r2 



In 



det^" 



-DIfAA) 



7r2/54 



det [-D\{A^)] det [-D\{A2)] 

X) 

dt{t^ - if/^Re [gBF,± (/?, mt, 0, v)] 
00 ^ 

^-[TVA(e*™i)+Tr^(e*"-^)], 



n=l 



(3.14) 



where 



gBF,+ = TrBFg{f3, mt, 0, u) = ^ g{P, mt, 0, vij - V2,k), (3.15) 



3,k: 



gBF- = TrBFg{/3,mt,0,v + vr) = g{j3,mt,0,vij - V2,k + vr). (3.16) 



We minimize the effective potential in eq. (3.14) with respect to the vij and V2^k and 



obtain the expected results. With antiperiodic (-) boundary conditions on fermions the 
deconfined phase is always favoured: 

vi,j = V2j = 0, \/j = l,...,N. (3.17) 

When periodic (+) boundary conditions are applied to bifundamental fermions 

TT vr 
vij = -, V2,j = --, V i = 1, N, (3.18) 

or any rotation, of all the eigenvalues, such that vij — V2^k remains an odd multiple of 



TT, in agreement with |47, 42]. Here the Z2 exchange symmetry between the gauge field 
eigenvalue angles vij and V2^k is clearly broken. However, the effective potentials are the 
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same regardless of the fermion boundary conditions. Again we look at the form of the 
partition functions: 

QCD{BF, +) : In Zbf,+ ~ g{vi,j - V2,k), vij = 7r/2, V2^k = -7r/2 

QCD{BF, -) : In Zbf- ~ givij - V2,k + tt), vij = V2,k = 0. 

(3.19) 

This is shown explicitly for = 8 in Figure (Right) where comparison is made with 
QCD(Adj,+/-). The conclusion in terms of the effective potentials is that ^Vbf,± = ^Adj- 
[where we considered U{N) QCD(Adj) for the purpose of comparing with the orientifold 
equivalence. For U{2N) QCD(Adj) VAdj- = '2Vbf,±] but there is no equivalence with 
VAdj,+ (at least not at one loop). The factor of 1/2 results because of the additional 
gauge boson term in QCD(BF) and the fact that the number of Dirac flavours multiplies 
the fermion term, where the same number of Majorana flavours (one Dirac flavour is two 
Majorana) multiplies the fermion term in QCD(Adj). Of course these conclusions also hold 
for the SU{N) theories in the large A^ limit. 

4. Conclusions 

We can summarize the results and draw some conclusions concerning one-loop calculations 
of the phase diagrams of QCD(R) for fermions in the fundamental, antisymmetric, symmet- 
ric, and adjoint representations. We have minimized the one loop effective potential with 
respect to the Polyakov loop eigenvalues, for a range of m/3, on the topology x M^, for 
small S^. For ABC on fermions the deconfined phase is favoured for all mf3 in SU{N) gauge 
theories with all the above fermion representations. With PBC on fundamental fermions 
a phase in which the Polyakov loop eigenvalue angles are all close to vr is preferred. For 
A^ odd this results in C-symmetry breaking since the eigenvalue angles only approach vr, 
including an extra 0{1/N) contribution. For PBC on antisymmetric/symmetric repre- 
sentation fermions a C-breaking phase is always favoured where the Polyakov loop angles 
prefer to be ib7r/2 [including additional 0{1/N) contributions for A^ not a multiple of 4]. 
For PBC on adjoint representation fermions and Nf = 1 Majorana fermion flavour the 
deconfined phase is favoured for all m/3 > 0. However, for Nf > 2 the phase structure is 
quite varied and depends on A". 

Concerning orientifold planar equivalence, if C-symmetry is broken in QCD(AS/S) 
by applying periodic boundary conditions to fermions, then the equivalence still holds 
because the effect of transforming boundary conditions and of C-parity breaking cancel 
each other out, causing the partition function to be the same [for the U{N) theory, and 
with 0{1/N) corrections for the SU{N) theory] as when antiperiodic fermion boundary 
conditions are applied. However, if C-symmetry is broken in some other way, such that the 
effective potential is not minimized to the same value as for QCD(AS/S,-) without C-parity 
breaking, then the equivalence fails. 

Considering orbifold equivalence, the breaking of Z2 exchange symmetry when periodic 
boundary conditions are applied to bifundamental representation fermions results in the 
same partition function as when antiperiodic boundary conditions are applied. Therefore 
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the equivalence still holds because the effect of the changing fermion boundary conditions 
from antiperiodic to periodic adds a factor of vr to the gauge fields that cancels out the 
exchange symmetry breaking. 

From our results, and from the derivation of the partition function in Appendix 
several well-established group theory equivalences are clear. For example, QCD(S) is the 
same as QCD(Adj) in SU{2) because the symmetric representation is the same as the 
adjoint (see Appendix^. QCD(AS) is the same as QCD(F) in SU{3) if — > — /i because 
the antisymmetric and the antifundamental representations correspond. 

Considering the partition function we also have, in U{N) [or SU{N) with N = 4, 8, 12, 
...], Zas- the same as Zas,+, even though the Polyakov loop eigenvalues differ. Of course, 
this is also true for Zs^- and .^5,+. The same is true for Zbf~ and Zbf,+ which are identical 
considering the U{N) theories [and agree to 0{1/N) corrections considering the SU{N) 
theories] and correspond to Z^g/g^^. in the large- limit. Whether these equivalences hold 
non-perturbatively is not clear. 

The results in this paper were calculated the case of /x = 0, however it is straightforward 
to include a finite chemical potential using the formulas in Appendix^ The only restriction 
is that the integral representation is only valid for /x < m, but the series representation 
often converges rather quickly. 
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Figure 5: N ^ 2: Vad.j(+) for Nf = 2 
Majorana or Vsymm(+) for A^/ — 1 Dirac 
flavour 



Figure 6: Vadj{+) for = 3 and Nf = 2 





Figure 7: Vadj{+) for = 4 and Nf ^2 Figure 8: Vadj(+) for = 5 and Nf ^2 
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Figure 9: Vadj{+) for A^ = 6 and Nf = 2 Figure 10: Vadj{+) for A^ = 7 and Nf = 2 
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Figure 11: Vadj{+) for = 8 and Nf ^2 Figure 12: Vad.j{+) for iV = 9 and iV/ = 2 
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Figure 13: Accessibility of the confined phase in QCD(Adj,+) for Nf — 2 and N from 2 - 19. The 
curve for ifn(3)^^^^ is a lower bound on the value of transition out of the confined phase. However, 
for each N, the point at mf3 = {''n(3)^^-^ + 0.1 does not belong to the confined phase. 
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Figure 14: Vadj{+) for N = (Left) Nf = 2; (Right) Nf = 3. These results suggest that 
increasing Nf causes the confined phase to be preferred up to a larger {mj3)crit- 
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Figure 15: Location of the transition to the deconfined phase phase in QCD(Adj,+) for N j from 
2 through 5, for aU N. 
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Figure 16: {XX)adj(+) for N = 2,Nf = 2 Figure 17: {XX)adj{+) for N = 3,Nf = 2 
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Figure 18: {XX)adj{+) for N = 4,Nf = 2 Figure 19: {XX)adj{+) for N = 5,Nf = 2 
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Figure 20: {XX)adj{+) for N = 6, Nf = 2 Figure 21: (AA)a£)j(-k) for N = 7, Nf = 2 
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Figure 22: (AA)adj(+) for N = 8,Nf = 2 Figure 23: (AA) az5j(+) for N^9,Nf^2 
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Figure 24: (Left) Orientifold planar equivalence: To check at one-loop that the equivalence is 
vahd for QCD(Adj,-), QCD(AS,+), QCD(AS,-), we plot the effective potential in all three cases 
for the U{N) theories at TV 8. For QCD(Adj) Nf^ = 1; for QCD(AS), = 1. Comparison 
with QCD(Adj,-|-) indicates that this theory is not equivalent to any of the others in this limit. 
The effective potential for QCD(Adj,-) is minimized in the deconfined phase where the Polyakov 
loop angles are Vi — Vi. This is also true for QCD(AS,-). For QCD(AS,-|-) Veff is minimized 
in the C-brcaking phase where Vi = tt/2 Vi. Regardless, the effective potentials of QCD(AS,+/-) 
correspond exactly. (Right) Orbifold equivalence: Analogously with the orientifold result we find at 
one-loop that orbifold equivalence is valid for QCD(Adj,-), QCD(BF,4-), and QCD(BF,-), but not 
QCD(Adj,+). Wc plot the effective potentials VAdj{±) with Nf = 1, ^Vbf{±} with ^ 1 for the 
U{N) theories with — 8. In QCD(BF,-) the effective potential is minimized in the deconfined 
phase with vi^i — W2,i — Vi. For QCD(BF,-|-) Veff is minimized in the Z2 interchange broken 
phase where vi^t = 7r/2 and W2,i — ~7r/2, Vi. However, the effective potentials Vbf,± correspond 
exactly, and are identical to 2VAdj,-- 
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A. Derivation of \n Zqcd{r) 



For Nf quark flavours in representation R, with masses mj, at chemical potential /ij, and 
at inverse temperature /3 = 1/T, the QCD(R) partition function is given by the (Euclidean 
space) path integral 



ZQCDiR) W,^) = J P^e-^-(^) J vme-^''^^''-<'^-~^-^^''>, (A.1) 

where the gauge field Af^{x) = Tg^^(x), /i = 0,1,2,3, the a indices are a = 0,...,N'^ — 
1, and the are the generators of SU{N) in the representation R. ^ is a vector of 
fermion fields containing Nf anti-commuting 4-spinors ^pf in representation R. ip contains 
the corresponding antifermion fields. M is the fermion mass matrix where {M)ffi = rufSff 
and is the fermion chemical potential matrix where (A^) J J' = fj^fSff. Ipj^ = 7^!?^ where 
Dfj^ is the covariant derivative 

D^ = df, + A^. (A.2) 

The are D{R) x D{R) matrices in the colour space, taken in the representation R of 
the fermion fields, and 7^ are the Euclidean space Hermitian gamma matrices satisfying 
{lfi-,lv} = 2(5^1,. We work in Euclidean space throughout this paper with the metric 
Qfiu = ^^Jl.v■ Note that in eq. ( |A.2| ) we rescaled A to absorb the coupling g according to 
gA^ A^. 

Sym{A) is the pure Yang- Mills theory action which contains only the boson contribu- 
tion. It is given by 

Sym = ^ j (^A*^^/^^) ' (A.3) 

where is the field strength defined by 

F^, = [D^, D,] = df,A, - dyA^ + [^^, A,\ . (A.4) 
Because the functional integrals in eq. (A.l) have an exponential that is quadratic in 



the fermion fields V'/ they are of the gaussian form and exactly solvable. The result is 

Zqcd{R) f^) = J-DA det (^^ - 70^^ + M) e-^^*^(^) . (A.5) 

We need the eigenvalues of Ipj^ to evaluate the determinant for fermion fields ip in 
representation R. For free fields this is more easily done in momentum space using 



\nZ = pV3 I ^lnZ{p), (A.6) 



d^p 

where now in eq. ( |A.5| ) is evaluated by having it act on the fermion fields in momentum 
space. These have the form 
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where V3 is the 3-volume, and lo^ are the Matsubara frequencies for periodic (+) or an- 
tiperiodic (-) boundary conditions apphed to fermions, 



cjj" = 2n7r//3, 

' (A.8) 
= (2n + l)7r//3. 

The compactification of also leads to periodicity in the gauge field A. In this case 
the boundary conditions are periodic (+) and we expand Ap, in a Fourier series according 
to 



^^(T,x) = i^e-"-A^ K, x). 



/5 



(A.9) 



Using eq. ( |A.2D for and eq. ( A.7 ) for the eigenvectors '0 we can compute the 
eigenvalues i\ of jp^ in frequency-momentum space. The result is 



%\ = i7ou;± + i7 • p + 7m^m^^ (A-IO) 

where i\ is a constant times the identity matrix in colour (k = 1,...,N), flavour (/ = 
1, ...,Nf), and spinor space. A^ are elements of the Lie algebra of SU{N). Here we have 
chosen A^ to he N x N and transform as the representation R. Then the T"" are the 
generators of SU{N) in the fundamental representation. 

Ipj^ is anti-Hermitian = —Ip^) since the gamma matrices are Hermitian, ojn and 
p are real, and the Afj_ are anti-Hermitian. Therefore, the eigenvalues iX of are pure 
imaginary and A G Also, since the Dirac operator satisfies {75, -^/j} = 0, and since 
{7O) 75} = 0, the eigenvalues of Ipj^ — 70/^/ show up in pairs in frequency-momentum space. 
The eigenvalue equations are 



(^^-7oA1)'0 = (iA-7oA1)V, ^^^^^ 
(^^ - 7o7W) (75V') = -iiX- loM) (75V') , 

which agrees with ||50| , |5l| for the fundamental representation with = 0. 

The fermion determinant in eq. (A. 5) is over the spinor indices, colour indices, flavour 



indices, as well as over frequency-momentum space. We first take the determinant over the 
colour indices {k = 1, ■■■,N) using the identity, Indet A = Trln^, for an aribtrary N x N 
matrix A. Then we take the determinant over the spinor indices (4 degrees of freedom for 
Dirac spinors and 2 for Majorana spinors), then frequency-momentum space indices, then 
finally over flavour indices. This proceeds as 
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lndet(^^-7oA^ + M) 

= Ttr In det (^^ - jqM + M) 

= 2 TYr ^ ^ In det - -foM + M) {-iX + -/qM + M)] 



n p 



2TVR^^lndet [- {iX - -foM)"^ + M^' 

n p 

2 TYi? ^ ^ In det [- (iw^ + Aq - M)^ - {ip + Af + M"^ 

n p 

r ^4 



(A.12) 



In 



[lo^ -iAo + i^f) + (p - iA) +m' 



To get from line 4 to line 5 we used, in A, the Euclidean gamma matrices 7^ in the chiral 
representation where 75 = 71727370 = diag{l, 1, —1, —1}. To see clearly the effect of 
antiparticles we can rearrange the result as follows 



I In [(u^ - Mo + i/x/)^ + (p - iAf + mj 



J^(ln Uutf + {u;f - if,f - Ao)y 



+ In 



(A.13) 



where ujf = y (p — iA) + mj. The second term in ( [A.13 ) is the antiparticle contribution 

which has ^0 ~^ ^0 and fif ^ —fJ-f with respect to the particle contribution as expected. 
To simplify later calculations we now take all the fermion masses and chemical potentials 
to be degenerate m/ = m {ujf = oj), fif = /u,V/. 

The final result for the fermion contribution is then 



In det {Ipj^ - 7o7W + M) 



(A.14) 



where u = fi — Aq. Defining the zero-temperature determinant to 1 and evaluating as in 
|3C] we get 



In det {lp^--foM + M) 



2V3Nf TrR 



d^p 
(2^ 



In (1 ^ e-^^^-^'A + In f 1 ^ e 



-I3{ui+u) 



(A.15) 



where the top sign corresponds to periodic boundary conditions [PBC (+)] on fermions 
and the bottom sign is for the usual case of antiperiodic boundary conditions [ABC (-)] 
on fermions. Spherical symmetry in p suggests evaluating the integral by converting to 
hyper-spherical coordinates using 
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and 



r(n + - 



rid/2) 



(2n)! 
n!22" 



TT. 



(A.16) 



(A.17) 



Then 



lndet(^^-7o7W + M) 



2V3Nf 



(27r)' 



i 



In f 1 T e-^('^-")') + In fl T e-^^'^+^^l 



(A.18) 



The Polyakov loop is defined as the path-ordered exponential of the temporal compo- 
nent of the gauge field, 

P{x) = Ve-^o d^Mx) _ (A. 19) 

For a constant background field defined by Ao = iv/P field the Polyakov loop is 



p = e>^^° = e^". 



(A.20) 



where we have chosen a gauge in which Aq is diagonal and v is real, diagonal and traceless 
with elements {v)ij = ViSij. Then 

P = diag{e^^S...,e^''^}. (A.21) 
We make a further simplification by taking = for z = 1, 2, 3. Then co = ^/■p^ -\-rrP- 

and 



lndet(^^-7oA4 -KM) 

= 2Vaf^f-^^ Tyr dp/-i [in (l T e-^('^-'*)p^) + In (l T e-^^'^+'^^P 

-i)/2R(d-i)/2 ^Rl^ „(d+i)/2 ^ +^ """P )K^d+i)/2{nPm). 



2id-3)/2^{d+l)/2pid-l)/2 



n=l 



(A.22) 

We are interested in the physical case of d = 3. Then using an integral representation 
for modified Bessel functions of the second kind: 

valid for v > —1/2 and Re{z) > 0, the log of the fermion determinant becomes 
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lndet(^^ -70-M + M) 



V^Nfm^ ^ (±1) 



n=l 



(A.24) 



./I 



dt(t2 _ 1)3/2 L ^ ^ ^t^^^ , 



where gR^± depends on the group representation of the fermions and the boundary condi- 
tions. We define the matrix 



giP, mt, ix,v) = J2 e-^"^-"^("^*-'^) 



n=l 



(A.25) 



gfS(mt-ti) _ 2 cos t; + e~/3(mt-^t) ' 
where the first hnc converges to the second hne for fi < ni. Wc study fermions in the 
fundamental (F), adjoint (Adj), symmetric (S), and antisymmetric (AS) representations. 
For periodic boundary conditions (+) gji^^ is defined according to gn^^ = Tr jig(P, mt, v). 
Then 



AT 



1=1 



AT 



9F,+ = Ttf g{l3, mt, fi,v) =Y g{(3, mt, n, Vi) = ^ 



-ivi _ p—l3{mt—iJ,) 



1=1 



e^("**-'*) - 2 cos Vi + e-^("**-'^) ' 



(A.26) 



N 



9Adj,+ = Tr^ 9{(i, rnt, n,v) = ^ g{l3, mt, ji, Vi - vj) - g{/3, mt, fi, 0) 



1 _ Q-P{mt-ii) 



(A.27) 



1^ f,l3{mt-ij) _ 2 cos {vi — Vj)+ e-Z^M-M) (=l3{mt-ij) _ 2 + e-l3{mt-ij) ' 



N 



gAS,+ = Ttas 9{I3, mt, fi,v) = ^ g{(3, mt, fi, Vi + vj) 

i<j=l 

^ p-i{vi+Vj) _ ^-f3(mt-j£) 



(A.28) 



E 

i<3= 



^ f,l3{mt-ij,) _ 2 cos(t;i + Vj)+ e-/^M-M) 



iV 

9S,+ = Tr5 g{P, mt, n,v) = ^ g{P, mt, /x, + vj) 

i<j=i 

^ Q-i{vi+Vj) _ p-l3{mt-fi) 

~ eP{mt-iJ,) _ 2cos(Di + Vj) + e-/'M-/') ' 



(A.29) 



where the form of the higher dimensional representations can be derived in terms of the 
fundamental using the Frobenius formula in combination with tensor product formulae from 
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Young tableaux as detailed in Appendix ^ To get the results for antiperiodic boundary 
conditions applied to fermions we just take v ^ v + n, so qr- = TrRg{l3,mt, iJ,,v + vr). 
Then 



N 



i=l 

N 

9Adj- = ^ g{P, mt, fi, Vi - Vj + vr) - g{(3, mt, fi, vr), 

N 

9AS-= 9{P,'mt,ij,Vi +Vj + tt), 

i<j=l 
N 

9s-= ^ g{P,mt,n,Vi + Vj + tt). 
i<j=i 



(A.30) 
(A.31) 
(A.32) 
(A.33) 



To obtain the phase diagram we need to minimize the free energy which is given by 
the effective potential: 



/ = Veff{m,fi) 



1 dlnZ 



-— In Z{m,fj,) 

1 



(A.34) 



[In ZB{m, n) + In Zpim, fx)] 



where the fermions have been integrated out and their contribution to the effective potential 
is: 



In Zp{m, fi) 



dt{t^ - 1)^/^ gR^± {P, mt, fi, v) +g]^^ {13, mt, -fi, v) 



(A.35) 



The one-loop contribution is calculated by introducing fluctuations yl^ around a slowly 
varying background field such that = + gA^. The free energy has the form 

f^T''[co + C2g^ + C3g^ + ...], (A.36) 
where the one-loop result corresponds to the cq term. The contributing Feynman diagrams 



at one loop are shown in Figure 25 



Considering just the one- loop contribution the gauge- fixed Lagrangian is 



^QCDiR) = [-5^,D\a)\ Al + [-^^(a)] + (^^ + M - ^^M) i>, (A.37) 
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i 




Figure 25: One-loop contributions to the free energy 



where C and C are the complex grassman ghost and antighost fields, respectively, which 
result from the gauge fixing. Then the path integral is 



QCDiR) {P,IJ) = j VAVCVC j PV'I?Ve"^o''d^/d'''^QCD{fl) 

V{gA)VCVC j V^Vi^e- ^0 I '^''''^QCD(R) ^ 



(A.38) 



So the problem is reduced to performing Gaussian integrals. The final result is 



Zqcd{r) = det [Ip^{a) +m- -fofif' det"^ [-D\d^{a)) . (A.39) 



The Yang-Mills theory result from the gluon and ghost contribution was calculated in 
|2S] and the result in terms of the Polyakov loop P = diagje**"!, ...,6^^^} is 



WM = -^lndet i-Dl^) 



oo ^ 



7r2/?4 ^ n' 



n=l 



(A.40) 



1 

pi 



1 ^ 

7-2 ^[Vi- Vj? (27r - [Vi - Vj]) 



247r2 



2 - (iV2 - 1) 

45 ^ ^ 



where [v] = v mod 27r. 

Now we can combine the boson and fermion contributions to get the full one-loop 
effective potential in terms of the Polyakov loop angles Vi. The result is 
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Vi-ioop{P,m,l3,tJ,) 
= -^lnZ{P,m,(3,ii) 



= ^ [-2Nf Indet (- {Dr{P) - ^if + rr?) + Indet {-D\^^{P)) 



m'^Nf ^ (±1) 



E 



n=l 



2 ^ 1 



n=l 



+ 



1 ^ 2 
^J:^- (2vr - - V,])' -^{N'- 1) 



If we take /i = then the result simphfies and we have 



(A.41) 



Vi-ioop{P,m,(3) 



2m^N 



7r2/32 



^ L_LRe [TrK(P«)] K2(n/?m) - ^ ^ -TV^l^'") 

n=l 



1 



dt{e - if/^Re [gR,± iP, mt, 0, v)] 



(A.42) 
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B. Higher dimensional representations 



In this appendix we show how to get higher dimensional representations of Tr/jP in terms 
of the fundamental TrpP where P E SU{N). In what follows we will refer often to Young 
tableau. A representation R is referred to in terms of its Young tableau by a comma 
separated list enclosed in parentheses {yi,y2, ■■■,yN-i)- Thus, an arbitrary representation 
R, with a Young tableau as shown in Figure 26, is represented by {yi,y2, ■■■,yN~i) = 
(p,g,r, ...). 



r boxes 



q boxes 



p boxes 



< N 
boxes 



Figure 26: Young tableau of an arbitrary representation R 

In what follows we will deal with representations that have either one or two rows of 
boxes in their Young tableaux, and so we will use the notation {p, q) for Young tableaux 
that have p columns with one row of boxes, and q columns with two rows of boxes. However, 
the procedure is valid in general for all representations. 

Higher dimensional representations R, with Young tableau {p, q), of SU (N) can be put 
in terms of the fundamental F, with Young tableau {N, 0), by using the Frobenius formula 



1 52, 53, In terms of the Polyakov loop the trace has the form 

Tr/?^ = ^,Y1 i^^^FPy' {TtpPy' ... (TrpPy- , (B.l) 

where n is the number of boxes in the Young tableau of the representation R, the sum is 
over all n! permutations j = {ii, j2i •••lin} of the symmetric group Sn, and xrH) is the 
group character, in the representation R, of the permutation j of the symmetric group Sn- 
These have a simple form when is a symmetric representation with Young tableau (n,0). 



This is given by [55| 



Tl 

nfc=i '^^"jk^- 

So for symmetric representations R with Young tableau (n, 0) 

Tr(n,o)^ = E n" I.. -I (T^fP)^'^ (Tr^P^)''^ ... {TrpPy- . (B.3) 
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After we have the symmetric representations it is simple to get the rest by evaluating 
tensor products using Young tableau. Alternatively, all the representations can be found 
from the Frobenius formula, but the characters in other representations are not as simple. 

Our procedure is best understood through example. Let us start with symmetric 
representations with Young tableau (n, 0). First we want to find all n! permutations j of 
the symmetric group Sn- We can obtain these from the elements of the symmetric group, 
Sn- The elements of Sn for n = 2, 3, 4 are 



S2 ■■ 


(1)(2) (12) 


S3: 


(1)(2)(3) (1)(23) 


(3)(12) 


(2)(13) 


(123) 


(132) 


54 : 


(1)(2)(3)(4) (1)(234) 


(1)(243) 


(2)(134) 


(2)(143) 


(3)(124) 




(3)(142) (4)(123) 


(4)(132) 


(12)(34) 


(13)(24) 


(14)(23) 




(1)(2)(34) (1)(3)(24) (1)(4)(23) (2)(3)(14) (2)(4)(13) (3)(4)(12) 




(1234) (1243) 


(1324) 


(1342) 


(1423) 


(1432) 



(B.4) 



Therefore the permutations j of Sn for n = 2, 3, 4 are 



n = 2 


{ii,j2} 


{2,0} {0,1} 


n = 3 


{ii,j2,i3} 


{3,0,0} {1,1,0} {0,0,1} 


n = 4 


{jl,j2,j3, ^4} 


{4,0,0,0} {1,0,1,0} {0,2,0,0} {2,1,0,0} {0,0,0,4} 



The permutations j can more simply be thought of as all possible solution vectors j = 
{jl,j2, -Jn} to 



lji+2j2 + -+nj„ = n (B.6) 

with > 0, Vi. 

At this point it is simple to get the higher dimensional representations Tr/jP from eq. 
(B.3) where i? is a symmetric representation (n,0). These are 



Tr(2,o)^' = I ^ 
Tr(3,o)^ = I 



{TrpPf + 3 (Tri^P) (Tr^^P^) _^ 2 (Tr^^P^) 



4! L 



(TVi^P)^ + 6 (Tr^P)^ {TtfP^) + 3 [TtpP'^f + 8 (Tr^P) (Tri^P^) + 6 (Tr^P^) 



(B.7) 

Now that we have a procedure for obtaining all the Tr^^P for symmetric representations, 
we can get the rest by means of tensor products. Some useful tensor product identities 
(which can be obtained from Young tableau) are 
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(1,...,O)0(1,...,O) = (2,...,0)©(0,1,...) 
(1,...,0)®(0,...,1) = (i,-,i)©(o,...,o) 
(2,...,0)55(0,...,1) = (2,...,l)e(l,...,0) 

(2, 0) (0, 2) = (2, 2) © (1, 1) © (0, 0) 
["..." indicates the appropriate number of zeroes.] 



(B. 



where we have brought back in the extra elements of {yi,y2, ■■■,yN-i) for the purposes of 
getting the adjoint (A) representation (1, 1). Also, other representations of interest are 
the symmetric (S) (2,...,0) and the antisymmetric (AS) (0,1,...). 

The first line of eq. (|B.8D allows us to get the antisymmetric representation from the 



symmetric, which we found in the first line of eq. 
From the first line of eq. (|B.8|) 



, TrsP 



Tr{o,i,...)^ 



(Tr(i,...,o)^') 



Tr(2,...,o)^ 



TtAsP = (TrFPy - TrsP 



(B.9) 



Tr^sP = 2 



iTrpPf - {TrpP^) 



The second line of eq. ( |B.8D allows us to get the adjoint representation from the 
fundamental (F) (1,...,0) and antifundamental (F) (0, ...,1) representations 



= (Tr(i,...,o)P) (TV(o,...,i)P) - 1 
Tr^P = |Tri.P|^ - 1. 



(B.IO) 



Other representations which are sometimes used are represented by (2, 1, ...) and (2, 2, , 
These are obtained from the third, and fourth line of eq. (p.8|), respectively. 



^H2,...,l)P 



(Tr(2,...,o)^') (Tr(o,...,i)P) - Tr^i,...,o)P 



TtfP^) liTipPf + (Tr i.p2 



TrpP, 



(B.ll) 



Tr(2,...,2)P 



[Tr(2,...,o)P) (Tr(o,...,2)P) 
^ ^TrpPf + {TipP'^) 



Tr(w)P-l 



TtfP^ + 



|TrirP|2. 



(B.12) 



These representations are often named by their dimension D which can easily be 
obtained using the "factors over hooks" rule D = F/H ||5^. For SU{5) and SU{4) the 
above representations have the following dimensions. 
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name 


(Vl, V2-, VN-l) 


o u yo J 


au (1)) 


fundamental (F) 




3 


4 


antifundamental (F) 


(0,...,1) 


3 


4 


adjoint (A) 


(1,...,1) 


8 


15 


symmetric (S) 


(2,...,0) 


6 


10 


antisymmetric (AS) 


(0,1,...) 


3 


6 


? 


(0,...,2) 


6 


10 


? 


(2,...,1) 


15 


36 


? 


(2,..., 2) 


27 


84 



For these representations the notation in the case of SU (2) needs to be considered specially 
since (yi, ...,yN-i) = (yi) is usually not sufficient to describe Young Tableaux for repre- 
sentations with Hi ^ where i > 1. Instead, wc can use (yi + y*) to count the columns of 
boxes in a single row where y^ is the number corresponding to yN-i above. This makes 
use of the fact that the fundamental and antifundamental representations are the same in 
SU{2). Hence, the adjoint representation is the same as the symmetric. Also, in SU{2) 
y2 represents only singlet contributions so the antisymmetric representation becomes the 
singlet. 
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